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The objective of this research was to study the effect of ecotype and planting location on thermal 
properties of big bluestem. Three big bluestem ecotypes (CKS, EKS, ILL) and a cultivar (KAW) were 
harvested in 2010 from four locations (Colby, Hays, and Manhattan, KS; and Carbondale, IL) and were 
evaluated for their specific heat, thermal conductivity, thermal stability, HHV (high heating value), and 
proximate contents. All populations revealed a large variation in specific heat (2.35—2.62 kJ/kg/K), 
thermal conductivity (77.85-99.06 x 10 -3 W/m/K), thermogravimetric analysis as weight loss during the 
heating process (71-73%), and HHV (17.64-18.67 MJ/kg). Specific heat of the big bluestem was signifi¬ 
cantly affected by planting location, ecotype, and interaction between location and ecotype. Planting 
location had stronger influence on specific heat than ecotype. Specific heat increased as temperature 
increased, and a linear correlation model for specific heat prediction was developed as a function of 
temperature. Ecotype, planting location, and the interaction of ecotype and planting location did not 
have a significant effect on thermal conductivity; however, density and particle size showed a completely 
opposite relationship on thermal conductivity. Both planting location and ecotype significantly affected 
HHV. Among all environmental factors, potential evapotranspiration had the most significant effect on 
thermal properties. 

Published by Elsevier Ltd. 


1. Introduction 

Renewable energy has received growing attention as people 
have become more conscious of the fossil fuel shortage and 
greenhouse gas emissions have been related to global warming [1 ]. 
In 2010, renewable energy resources supplied 8% of the nation’s 
total energy consumption and up to 8.05 quadrillion Btu (increased 
from 6% in 2009). Biomass led the other renewable energy re¬ 
sources (such as wind, solar, geothermal, and hydroelectric) by 
contributing 53% of the nation’s renewable energy supply in 2010 
[1]. Biomass resources include various natural and derived mate¬ 
rials, such as woody and herbaceous species, wood wastes, bagasse, 
agricultural and industrial residues, waste paper, municipal solid 
waste, sawdust, biosolids, grass, waste from food processing, 
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animal wastes, aquatic plants, algae, etc. [2]. Big bluestem ( Andro- 
pogon gerardii ) is regarded as second-generation biomass and 
recently has been proposed as promising energy crops because 
their growth requires few agricultural inputs (fertilizer and pesti¬ 
cides). Propheter and Staggenborg reported that only one-tenth of 
nitrogen filterlizers were required for big bluestem during the first 
two years comparing with annual biofuel crops such as corn, soy¬ 
bean and sorghum [3]. In addition, those annual biofuel crops had 
much higher soil nitrogen removal rate (~160 kg ha -1 ) than big 
bluestem (~30 kg ha -1 ) [3]. Lynd reported that perennial grasses 
such as big bluestem require substantially less fertilizer and pes¬ 
ticides than corn as biofuel feedstocks [4]. In the United States, big 
bluestem dominates the tallgrass prairie of North America and is a 
major component of prairie biomass [5,6]. Moreover, tolerance to 
heat and drought have enabled big bluestem to fill the deficiency in 
grasslands in the Midwestern U.S. when cool-season grasses (C3) 
are unproductive [7-9]. Recently research reported that succes- 
sional herbaceous vegetation, such as big bluestem and alfalfa, on 
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marginal lands in the Midwestern U.S. states can not only provide 
greenhouse gas emissions mitigation, but also produce substantial 
proportion of future biofuel energy target [10]. 

Thermal, biological, and physical processes are the three major 
technologies that help make use of a wide variety of biomass. In 
thermal conversion technologies, direct combustion and co-firing 
with coal were first utilized for electricity production and once 
were responsible for over 97% of the world’s bioenergy production 
[11 ]. Pyrolysis has attracted the highest interest because it produces 
bio-oil, which can be used as a fuel in stable engines and converted 
into chemicals such as bio-lime nitrogen fertilizer [12]. Biomass 
gasification has been researched extensively due to its higher effi¬ 
ciencies compared with combustion, and fast pyrolysis is still at a 
relatively early stage of development [13]. Torrefaction, another 
promising thermal process, improves the quality of biomass in 
terms of heat content, physical properties, and chemical composi¬ 
tion for combustion and gasification applications [14]. 

Understanding, predicting, and controlling these thermal pro¬ 
cesses and designing processing equipment require knowledge of 
the thermal properties of biomass, such as specific heat, thermal 
conductivity, thermal stability, HHV (high heating value), and 
proximate contents. 

Specific heat of a substance (kJ/kg/K) is defined as the amount of 
heat required to increase the temperature of a unit of mass by one 
degree. Specific heat affects the total energy required for thermal 
conversion of biomass into biofuels. The specific heat of biomasses 
depends largely on their composition; using the specific heat of 
each component of a mixture and the mass fraction is usually 
sufficient to predict the specific heat of the mixture [15]. Although 
the method based on the specific heat of the components in the 
mixture is most widely used to predict specific heat because of its 
simplicity, experimentally determined value is usually higher than 
predicted value [16]. Koch utilized DSC (differential scanning 
calorimetry) as a convenient technique for measuring the specific 
heat of wood and bark of 72 spruce pine trees [17]. 

Thermal conductivity of a material (W/m K) is a measure of its 
ability to conduct heat. Thermal conductivity of biomass depends 
mostly on composition and the characteristics of the biomass that 
affect the heat flow paths through the material. Mohsenin reviewed 
thermal conductivity measurement techniques for both steady- 
state and transient-sate transfers [18]. The heated probe method 
is simple, fast, and requires only a small sample, and it has been 
widely used for thermal conductivity determination [19]. Thermal 
stability, the ability of a material to resist changes in physical shape 
or size as its temperature changes, is essential to understand and 
predict the reactions and kinetics during biomass thermal conver¬ 
sion. TGA (thermogravimetric analysis) is the usual technique for 
determining thermal stability by quantitative measurement of 
weight changes (loss/gain) associated with thermally induced 
transition as a function of temperature or time [20]. High heating 
value (MJ/kg) is an important thermal parameter to characterize 
the amount of energy produced by the combustion of a unit mass of 
a material. Proximate analysis is a simple and rapid procedure for 
defining the substance energy content and determining how clean 
and efficient the substance is for the purpose. 

Our recent research on big bluestem showed that planting 
location and ecotype as well as interaction of planting location and 
ecotype had significant effects on chemical and elemental 
composition of big bluestem [21]. In addition, our study also 
showed that bio-oil yield from big bluestem through hydrothermal 
conversion was significantly affected by both ecotype and planting 
location [22]. However, we found no research on the thermal 
properties of big bluestem, especially the effects of ecotype and 
planting location on the thermal properties of big bluestem. 
Therefore, the objectives of this research were to characterize the 


thermal properties of big bluestem and to study the effects of 
ecotype and planting location on thermal properties of big blue¬ 
stem and thus, fill a critical gap in fundamental knowledge of 
thermal properties of a valuable bioenergy grass. 

2. Materials and methods 

2 A. Materials 

Three big bluestem ecotypes, CKS (Cedar Bluffs [CDB]), EKS 
(Konza [ICON]), and ILL (12Mile [12M]), and the KAW cultivar, which 
is widely planted to restore marginal lands, were harvested from 
reciprocal garden plots in four planting locations (Colby, Hays, and 
Manhattan, KS; and Carbondale, IL) in 2010. Among the four loca¬ 
tions, the Colby planting site was used to test the threshold of 
drought tolerance and the possibility for planting in drier Great 
Plains locations. Two populations from each ecotype were evalu¬ 
ated for thermal properties. No nitrogen fertilizers and pesticides 
were applied to the big bluestem. Details of seed collection and 
planting location have been described previously [21 . The big 
bluestem samples were ground into powder using a RetschSM2000 
cutting mill (Haan, Germany) with a 1.0-, 2.0-, and 4.0-mm sieve, 
respectively. For thermal conductivity measurement, only 2.0- and 
4.0-mm particle sizes were used. After grinding, each sample was 
fully mixed in sealed plastic storage bags. To eliminate any error 
that might be caused by water, samples were dried at 45 °C for 24 h 
and saved in the plastic bag before measuring specific heat, thermal 
conductivity, high heating value, and proximate content. All mea¬ 
surements were carried out in short time at room temperature of 
22 °C and 40% humidity. 

2.2. Specific heat by differential scanning calorimetry (DSC) 

Specific heat of big bluestem was measured with DSC Q200 
V24.4 instrument (TA Instruments, New Castle, DE) calibrated with 
indium and zinc. An empty sealed pan was used as a reference for 
every measurement. Three-step scans were carried out in this 
study. The first scan was conducted with an empty hermetic pan to 
determine the baseline background heat flow, which was sub¬ 
tracted from subsequent measurements. Next, sapphire was 
weighed and sealed in a pan for determination value of E. E was the 
calibration constant and calculated by using Equation 1 : 

Cp S x H r x M 

E= H x 60 (1) 

where Cps is specific heat of sapphire, which was standard and 
obtained from the literature (kJ/kg/K); H r is heating rate, which was 
10 in this study (K/min); M is sapphire mass (mg); H is measured 
heating value (mW); and 60 is conversion constant (minute to 
second). 

For the sample run, an empty pan and a pan with a 5-mg sample 
were placed into the DSC. The specific heat of the sample is 
calculated by transposition equation and substituting E. Large- 
volume stainless steel pans were used. All measurements were 
held at 323 K for 10 min, scanned from 323 K to 473 K at a heating 
rate of 10 K/min, and then held at 473 I< for another 10 min. The 
sample was characterized in an inert environment by using nitro¬ 
gen with a gas flow rate of 50 ml/min. 

2.3. Thermal conductivity by probe method 

Fig. 1 shows a diagram of the experimental apparatus for mea¬ 
surement of thermal conductivity using the heated needle probe. 
The container is filled with biomass with a confined particle size 
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Fig. 1. The diagram of apparatus for thermal conductivity measurement using the 
probe method. 


2.5. High heating value 

Gross energy contents were determined by means of a calo¬ 
rimeter (IKA-Calorimeter C 200, IKA-Werke GmbH and Co. KG, 
Staufen, Germany) with a benzoic acid standard. About 1.00 g of 
each pelleted sample was put into an adiabatic bomb and burned to 
ash. All samples were ground by a miller with a 1.0-mm sieve. 
Powder samples were compacted into pellets for measurement to 
reduce error caused by incomplete combustion resulting from dry, 
loose samples blown away during sudden release of volatiles. In 
addition, heating value was calculated based on the elemental 
composition of the biomass sampled. 

2.6. Proximate analysis 

Ash content was determined according to ASTM D-1102-84, 
“Standard test method for ash in wood” [23]. Volatile matter 
determination was made in accordance to ASTM E872-82, “Stan¬ 
dard test method for volatile matter in the analysis of particulate 
wood fuels” [24]. Fixed carbon was determined by subtracting the 
summation of the moisture, volatile matter and ash contents from 
the total sample mass. 


and sample density. The straight needle probe (60 mm; Thermal 
logic, Washington, USA) containing a heater wire as heat source 
and a thermocouple as a temperature-measuring device is inserted 
at the center of the container. The container’s diameter is 30 mm to 
ensure measurement time (3 min) shorter than the time allowing 
heat transfer to reach the wall of the container. After the initial 
temperature is equilibrated at room temperature, the heating wire 
is activated and heated at a constant rate of energy input supplied 
by a pair of AAA batteries. Then, the temperature rise overtime is 
measured with thermocouple and analyzed by Data Acquisition 
System NI USB-9161 (National Instruments). For the mathematical 
analysis, thermal conductivity can be expressed by Equation 2: 


k = 


Mi) 

4iz(T - T 0 ) 


( 2 ) 


2.7. Statistical analysis 

The reported thermal property values of big bluestem samples 
are the average of at least two replicates. Data were analyzed with 
ANOVA (analysis of variance) and Tukey’s studentized range 
(Honestly Significant Difference HSD) test in SAS (SAS Institute, Inc., 
Cary, NC). In general, fully balanced ANOVA tests were performed 
by following the GLM (general linear model) procedure. Correla¬ 
tions were determined using stepwise multiple regressions and 
multiple linear regression analysis. In addition to investigating ef¬ 
fects of environmental conditions on thermal properties of big 
bluestem, we performed multiple univariate linear regression 
analysis. For all weather data, we used the National Oceanic and 
Atmospheric Administration database and historical records dating 
back to 1961 (see [21] for further details). 

3. Results and discussion 


where k and q are the thermal conductivity (W/m/K) and heater 
power dissipated per unit length (W/m), respectively. T 0 and to are 
temperature (I<) and time (s) at initial condition. T and t represent 
temperature (I<) and time (s) since the probe heated is energized. 
The heat power q is calculated using Equation 3: 


3.1. Specific heat 

Fig. 2 shows the typical DSC plots of big bluestem, in which heat 
flow was investigated as a function of temperature. The maximum 


q = SqKm) (3) 

^ref 

where Rm is the known probe heater resistance (1041.5 Q) and F re f 
(1 Q) and ft re f (0.04 V) are voltage and resistance in the reference 
resistor, respectively. 

2.4. Thermogravimetric analysis 

The determinations were performed with the Perkin—Elmer 
Pyrisl thermogravimetric analysis (TGA) instrument (Norwalk, CT) 
to record the sample mass change with temperature over the 
course of the paralysis reaction. Initial TGA measurement was 
carried out at a heating rate of 10 °C/min up to 900 °C to establish 
the temperature range required for the following investigations. 
About 8 mg of each sample obtained from a tensile bar was placed 
in the pan and heated from 30 to 600 °C at a heating rate of 30 °C/ 
min in a nitrogen environment. 


- Baseline0606.001 

- TOW-111.001 

- RF0606.001 
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Fig. 2. Typical calorimeter plots for specific heat determination of big bluestem (green 
curve shows the baseline, and red and blue curves represent the DSC curve of big 
bluestem and a reference). 
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Fig. 3. Specific heat of big biuestem as a function of temperature. 

temperature of DSC measurement was set at 473 K to prevent ex¬ 
plosion of the sealed pan from emissivity of sample pyrolysis. A 
previous study did not recommend DSC analysis without the lid, 
although this approach can heat to higher temperature and was 
more approximate for biomass pyrolysis and gasification [25]. The 
green curve is the baseline with the sealed pan at a heating rate of 
10 K/min and was subtracted from following measurements. The 
red curve and blue curve represent the biomass sample and refer¬ 
ence (sapphire) under the same conditions, respectively. The heat 
flow curve includes the heats evolved during decomposition and 
reaction of the pyrolysis sample material. Therefore, during an 
endothermic heat effect, specific heat values increase, whereas 
specific heat decreases during an exothermic reaction [26]. The 
heat flow curve of the sample shows a decreased trend after 420 K, 
although three curves exhibit the positive relationships between 
heat flow and temperature from 323 I< to 420 K. Because the 
samples were previously dried to eliminate error from water, this 
endothermic reaction was associated with decomposition of 
hemicellulose [27]. Cultivar KAW planted in IL exhibited the most 
pronounced endothermic peak in this range due to its relatively 
higher specific heat value and hemicellulose content. The typical 
relationship between the specific heat and temperature of big 
biuestem is illustrated in Fig. 3. The value of specific heat increased 
linearly from 1.73 to 2.24 kJ/kg/K as temperature increased from 


323 K to 420 I< and followed a second-order polynomial in the 
whole range. This result suggests that the specific heat of big 
biuestem increased with heating until the component decom- 
postion. The average value of the specific heat of big biuestem at 
420 K is 2.46 kJ/kg/K (Table 1), which is in line with the previously 
reported range of 2.30-2.60 kJ/kg/K for wood and grass particle 
biomasses 28]. Stepwise multiple regression analysis was con¬ 
ducted to relate specific heat to temperature with factors I< and I< 2 . 
The best equation accounted for 85.1% of variation with a standard 
error of the estimate of 0.012 expressed as Equation 4: 

Specific heat of big biuestem = -0.1581T 2 +124.64T-22281 

(4) 

Equation 4 agreed well with the work of Gupta et al. [29 , which 
indicated that the specific heat of biomass increased by 40% and 
basically followed a linear pattern before 430 I< then followed a 
second-order polynomial for the entire pyrolysis process (Fig. 3). 

The bars on left side of Fig. 4 show the effects of planting loca¬ 
tion on the specific heat of big biuestem. Big biuestem populations 
planted in Illinois generally had a higher specific heat at 420 K, with 
an average of 2.62 kJ/kg/K compared with averages of 2.34 kJ/kg/K, 
2.40 kJ/kg/K, and 2.35 kJ/kg/K for populations planted in Colby, 
Hays, and Manhattan, KS, respectively. The average specific heat of 
big biuestem planted in Illinois was 0.25 kJ/kg/K, which is higher 
than those from Kansas and indicates that the same big biuestem 
populations would have around 10.4% higher value of specific heat 
if planted in Illinois instead of Kansas. The effect of ecotype on 
specific heat is shown on right side of Fig. 4. CKS, EKS, and ILL 
ecotypes had statistically similar average specific heat values 
regardless of planting location, but native cultivar KAW had 
significantly higher specific heat value than other ecotypes, as 
indicated by the different letters analyzed by Tukey’s HSD test. This 
result was in agreement with previous studies [21,22]. Table 2 
shows that the interaction between planting location and ecotype 
also had significant effect on specific heat of big biuestem 
(p < 0.05); however, the interaction effect was much less than that 
of planting location and ecotype and had a smaller F- value (11.3). It 
is noted that the role of the location was always a greater source of 
variation than ecotype and the interaction between ecotype and 
location for specific heat value based on larger F -value (Table 2). 

3.2. Thermal conductivity 

Table 3 summarizes the calculated thermal conductivity of big 
biuestem at three different densities and particle sizes at 298 K. As 


Table 1 

Specific heat, thermal conductivity, and thermogravimetric parameters by planting location and ecotype. 


Sample 

Specific heat (kJ/kg/K) 

Thermal conductivity 
(mW/m/K) 

Lonset (°C) 

Tend (°C) 

Tmax (°C) 

Weight loss (%) 

CKS-Colby 

2.23 ± 0.05 

79.28 ± 2.02 

296.1 ± 1.2 

404.9 ± 7.2 

357.6 ± 7.5 

71.53 ± 0.1 

EKS-Colby 

2.33 ± 0.05 

80.49 ± 2.83 

291.8 ± 2.6 

410.3 ± 2.8 

363.9 ±1.7 

71.35 ± 0.1 

ILL-Colby 

2.42 ± 0.05 

90.43 ± 1.77 

296.2 ± 7.3 

406.3 ± 0.7 

356.2 ± 2.1 

71.70 ± 0.4 

KAW-Colby 

2.50 ± 0.06 

87.74 ± 1.94 

301.2 ± 0.2 

412.1 ± .02 

364.1 ± 0.1 

74.33 ± 0.5 

CKS-Hays 

2.29 ± 0.06 

89.02 ± 3.96 

292.4 ± 0.8 

405.9 ± 1.8 

363.0 ± 1.0 

71.13 ± 0.6 

EKS-Hays 

2.32 ± 0.05 

82.00 ± 2.38 

291.7 ± 1.1 

407.6 ± 7.7 

356.6 ± 0.5 

71.25 ± 2.7 

ILL-Hays 

2.42 ± 0.02 

85.48 ± 1.15 

291.5 ± 5.4 

410.2 ±4.4 

364.6 ± 5.3 

72.08 ± 2.2 

KAW-Hays 

2.54 ± 0.02 

82.96 ± 2.23 

290.4 ± 0.5 

407.3 ± 0.4 

364.1 ± 0.1 

72.38 ± 0.5 

CKS-MH 

2.43 ± 0.04 

85.49 ± 3.55 

299.2 ± 1.7 

404.1 ± 1.7 

353.2 ± 0.7 

70.95 ± 0.2 

EKS-MH 

2.38 ± 0.03 

84.05 ± 0.80 

302.9 ± 2.9 

407.3 ± 3.5 

358.4 ± 2.5 

70.18 ± 2.0 

ILL-MH 

2.15 ± 0.03 

83.31 ± 3.07 

299.7 ± 4.8 

409.5 ± 4.5 

366.2 ± 0.7 

72.95 ± 1.1 

KAW-MH 

2.44 ± 0.06 

90.04 ± 3.78 

304.1 ± 0.1 

411.4 ±0.5 

362.1 ± 0.1 

72.28 ± 0.4 

CKS-ILL 

2.44 ± 0.06 

82.22 ± 3.071 

299.5 ± 2.5 

401.2 ± 2.8 

354.9 ± 0.6 

73.63 ± 0.3 

EKS-ILL 

2.56 ± 0.04 

87.76 ±3.11 

300.9 ± 4.4 

399.0 ± 3.3 

352.8 ± 3.8 

73.28 ± 1.1 

ILL-ILL 

2.62 ± 0.04 

87.58 ± 0.04 

299.5 ± 5.8 

394.6 ± 0.1 

353.4 ± 0.8 

72.63 ± 0.3 

KAW-ILL 

2.85 ± 0.04 

96.39 ± 3.78 

298.5 ± 0.7 

405.1 ± 0.1 

359.3 ± 0.4 

74.38 ± 0.5 
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Fig. 4. Comparisons of average specific heat of big bluestem as affected by planting 
location and ecotype. Different letters (a and b) above the standard deviation bars 
indicate they are statistically different at a 95% confidence level. 


density increased from 300 to 360 kg/m 3 , the value of thermal 
conductivity of big bluestem with particle sizes of 1 mm, 2 mm, and 
3 mm increased by 32%, 44%, and 49%, respectively. This result is 
probably explained by the fact that compression decreases the 
proportion of air held in samples, and air has lower thermal con¬ 
ductivity (24 x 1CT 3 W/m/K) than biomass (~85 x 1CT 3 W/m/K). 
The dependence of biomass thermal conductivity on sample den¬ 
sity was also reported by previous research [30]. Moreover, a strong 
negative correlation between thermal conductivity and particle 
size was reported by Hankalin et al. [31]. Biomass with smaller 
particle size has a higher thermal conductivity, which results in a 
better heat transfer and/or mass transfer in thermal processing. 
Thus, biomass with small particle size is considered advantageous 
to thermal conduct efficiency and has a higher production yield in 
the thermal conversion process [32]. 

To study the effects of planting location and ecotype on thermal 
conductivity of big bluestem, the probe method was conducted 
with the 1-mm sample in 360 kg/m 3 at 298 K. Thermal conductivity 
values were in the range of 77.85-99.06 x 10 -3 W/m/K at room 
temperature, depending on planting location and ecotype (Table 1 ). 
This result was higher than the thermal conductivity of woody and 
grass biomass (46.12-76.23 x 10 -3 W/m/K) reported by Yang et al. 


Table 2 

Effects of location, ecotype, and interaction between planting location and ecotype 
on the thermal properties of big bluestem. 


Thermal properties 

Source of Location 
variation 

Ecotype 

Location x ecotype 

Specific heat 

F 

65.87 

45.63 

11.29 


P 

<0.01 

<0.01 

<0.01 

Thermal conductivity 

F 

1.00 

0.35 

0.94 


P 

0.42 

0.79 

0.52 

Lonset 

F 

13.11 

0.52 

1.09 


P 

<0.01 

0.67 

0.42 

lend 

F 

10.64 

2.92 

1.38 


P 

<0.01 

0.07 

0.27 

Tnax 

F 

10.05 

6.10 

4.43 


P 

0.01 

0.01 

0.01 

Weight loss 

F 

4.71 

4.07 

1.00 


P 

0.02 

0.03 

0.48 

High heating value 

F 

43.02 

11.43 

1.26 


P 

<0.01 

<0.01 

0.33 

Volatile matter 

F 

313.13 

182.43 

33.68 


P 

<0.01 

<0.01 

<0.01 

Fixed carbon 

F 

50.16 

5.51 

1.53 


P 

<0.01 

0.01 

0.22 

Ash 

F 

1426.49 

278.24 

85.55 


P 

<0.01 

<0.01 

<0.01 


Table 3 

Effects of density and particle size on thermal conductivity of big bluestem. 


Density (kg/m 3 ) 

Thermal conductivity (mW/m/K x 

10" 3 ) 

Particle size (mm) 


1.0 

2.0 

3.0 

300 

63.58 

57.89 

53.69 

330 

79.90 

75.95 

72.74 

360 

84.08 

83.21 

80.04 


[27] and is partly explained by the fact that the thermal conduc¬ 
tivity was measured on samples with higher sample density in this 
study. In Fig. 5, although there was no significant difference in 
mean value of thermal conductivity between the planting locations 
and ecotypes, the thermal conductivity of the Illinois location 
(average 88.49 x 10“ 3 W/m/K) and KAW (average 86.70 x 10“ 3 W/ 
m/I<) was higher than those of the Colby location (average 
74.999 x 10 3 W/m/K) and CKS ecotype (average 78.55 x 10“ 3 W/ 
m/K). It is suggested that location and ecotype did not have sig¬ 
nificant effects on the average thermal conductivity from west to 
east, probably because thermal conductivity was not affected solely 
by location and ecotype. Table 5 showed the high correlation 
( R 2 = 0.99, 0.98, and 0.95) between potential evapotranspiration, 
aridity index, precipitation since 1961, and thermal conductivity, 
indicating that these three environmental factors may play a sig¬ 
nificant role in thermal conductivity values of big bluestem. 

3.3. Thermogravimetric analysis 

The thermogravimetric data were summarized in Table 1. Fig. 6 
shows the typical TGA and DTG (Derivative Thermogravimetric 
Analysis) curves of big bluestem. In general, all big bluestem sam¬ 
ples revealed similar curves: samples lost 10% of their weight 
within the temperature range of 30-250 °C, and significant weight 
loss occurred between 250 °C and 400 °C. The weight loss of big 
bluestems growing in Kansas and those from Illinois were 71% and 
73%, respectively, after achieving thermogram plateaus during 
heating. This result suggested that big bluestem growing in Illinois 
contained a higher amount of thermal decomposition contents 
(cellulose, hemicellulose, and lignin) compared with samples from 
the Kansas growing locations. The weight loss of big bluestem was 
significantly affected by both planting location and ecotype 
(Table 2).This result followed a similar trend in previous chemical 
composition analysis. In other words, big bluestem from Illinois and 



Fig. 5. Comparisons of average thermal conductivity of big bluestem as affected by 
planting location and ecotype. Different letters (a and b) above the standard deviation 
bars indicate they are statistically different at a 95% confidence level. 































































Table 4 

High heating value and proximate analysis of big bluestem. 
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HHV (MJ/kg) 

Prediction l a (MJ/kg) 

Prediction 2 b (MJ/kg) 

Volatile matter (%) 

Fixed carbon (%) 

Ash (%) 

CKS-Colby 

18.02 ± 7.36 

17.92 

17.73 

74.27 ± 0.2 

17.57 ± 0.2 

8.16 ± 

EKS-Colby 

17.45 db 40.35 

17.28 

17.70 

73.07 ± 0.2 

17.51 ± 0.1 

9.41 ± 

ILL-Colby 

17.28 db 7.59 

16.29 

17.54 

74.84 ± 0.2 

17.40 ± 0.1 

7.76 ± 

KAW-Colby 

18.39 db 23.11 

18.14 

18.43 

76.72 ± 0.3 

17.02 ± 0.4 

6.26 ± 

CKS-Hays 

17.90 db 6.52 

18.26 

17.51 

73.29 ± 0.3 

17.94 ± 0.4 

8.77 ± 

EKS-Hays 

17.60 ±21.30 

17.86 

17.20 

71.20 ±0.1 

17.50 ± 0.3 

11.31 ± 

ILL-Hays 

17.98 ± 28.40 

17.47 

17.59 

73.40 ± 0.3 

17.58 ± 0.3 

9.02 ± 

KAW-Hays 

17.96 ± 2.87 

18.92 

18.63 

74.85 ± 0.3 

17.54 ± 0.2 

7.61 ± 

CKS-MH 

18.29 ±37.63 

19.08 

17.76 

73.50 ± 0.4 

17.41 ± 0.7 

9.10 ± 

EKS-MH 

18.02 ± 101.65 

18.92 

17.59 

73.41 ± 0.1 

17.59 ± 0.3 

9.01 ± 

ILL-MH 

18.38 ± 23.17 

19.61 

17.92 

75.55 ± 0.2 

17.53 ± 0.1 

6.92 ± 

KAW-MH 

18.29 ±39.00 

18.92 

18.56 

74.48 ± 0.1 

17.28 ± 0.1 

8.24 ± 

CKS-ILL 

18.57 ± 45.56 

20.14 

17.93 

75.67 ± 0.1 

18.88 ± 0.1 

5.45 ± 

EKS-ILL 

18.51 ± 17.16 

19.69 

17.81 

76.48 ± 0.1 

19.06 ± 0.2 

4.46 ± 

ILL-ILL 

18.61 ± 2.18 

19.29 

17.66 

75.75 ± 0.1 

19.21 ± 0.1 

5.04 ± 

KAW-ILL 

19.05 ± 6.14 

20.16 

18.45 

77.51 ± 0.13 

18.10 ± 0.1 

4.39 ± 

RSS C 

NA 

11.4 

4.0 

NA 

NA 

NA 


a HHV = 339 C + 1214(H—0/8) 226 H + 105 S [36]. 
b HHV = 353.6FC + 155.9VM - 7.8ASH [37]. 

c RSS(residual SUm Of squares) = E?=l ( HHV predicted ~ HHV me asurement) 2 - 


KAW ecotypes had higher carbohydrate and lignin contents with 
weight loss than their counterparts. 

For DTG profiles, the rate of weight loss versus temperature of 
big bluestem, all curves exhibited three-stage thermal decompo¬ 
sition behaviors. The first stage occurred at peak temperature 
around 80-85 °C, which was related to the evaporation of water in 
the samples prior to the bulk of the weight loss. The second peak, 
decomposition of hemicellulose, was relatively small and located at 
195-255 °C, which is partly overlapped by the biggest peak. The 
third stage appeared from 200 to 510 °C, with 65% weight loss and 
peak temperature around 355—365 °C, indicating thermal decom¬ 
position of the main constituents of big bluestem, including cellu¬ 
lose and lignin. Hemicellulose was easy to decompose at a low 
temperature of 200-315 °C because of its random, amorphous, and 
rich branch structure [27]. The thermal decomposition of cellulose 
occurred at temperatures ranging from 275 to 350 °C 33]. Among 
the three constituents, lignin was the most difficult to decompose 
under the whole temperature ranging from 200 to 600 °C because 
its complicated oxygen functional groups have different thermal 
stabilities [34 . Based on statistical results in Table 2, although no 
significant difference in averages of T one and T en d was observed 
among ecotype, those thermogravimetric parameters of big blue¬ 
stem were significantly affected by planting location (p < 0.05). 
T max was significantly affected by all three factors. These results 
suggested that planting location rather than ecotype and interac¬ 
tion influenced the thermal stability of big bluestem. Furthermore, 
the weight loss of big bluestem gradually increased to about 3% as 


planting location changed from west to east. This result qualita¬ 
tively confirms the variation trend of contents of carbohydrates and 
lignin in previous research 21 ]. Multiple univariate linear regres¬ 
sion analysis was conducted to find significant associations of 
environmental predictors with thermogravimetric parameters, and 
the results showed 89% of the variation in the end time of thermal 
decomposition (T en d), 86% of the variation in the maximum time of 
thermal decomposition (T max ), and 71% of the variation in the 
weight loss were explained by potential evapotranspiration. Pre¬ 
cipitation in 2010 also explained a large variation in the onset time 
of thermal decomposition (T on set) of the big bluestem samples with 
coefficients of determination ( R 2 = 0.70) in growing year 2010 
(Table 5). 

3.4. High heating value and proximate analysis 

Heating value is the most important parameter in characterizing 
a substance as combustible, and it is widely used in the determi¬ 
nation of a number of additional thermal properties such as 
enthalpy of formation and adiabatic reaction temperature. High 
heating value (HHV), representing the heat of combustion relative 
to liquid water as the product, was measured with an adiabatic 
oxygen bomb calorimeter and compared with the predicted HHV 
based on ultimate and proximate analysis methods (Table 4). The 
proximate and ultimate analyses of biomass are essential for their 
efficient and clean utilization, whereas the HHVs of these materials 
determine the quantitative energy content of these fuels. The HHV 


Table 5 

Effects of environmental conditions on thermal properties of big bluestem analyzed by the multiple univariate linear regression analysis. 


Thermal properties 

PPT a 2010 (cm) 

PPT a since 1961 (cm) 

GDD b since 1961 (cm) 

GDD b 2010 (cm) 

PET C (cm) 

Aridity index 

Specific heat 

0.21 

0.62 

0.14 

0.44 

0.80 

0.70 

Thermal conductivity 

0.61 

0.95 

0.47 

0.63 

0.99 

0.98 

Tmset 

0.70 

0.56 

0.30 

0.07 

0.42 

0.49 

Fend 

0.32 

0.75 

0.19 

0.43 

0.89 

0.81 

Tnax 

0.38 

0.77 

0.15 

0.25 

0.86 

0.80 

Weight loss 

0.14 

0.54 

0.05 

0.23 

0.71 

0.61 

High heating value 

0.69 

0.98 

0.50 

0.59 

0.99 

0.99 

Volatile matter 

0.22 

0.56 

0.03 

0.07 

0.64 

0.58 

Fixed carbon 

0.28 

0.70 

0.23 

0.55 

0.86 

0.77 

Ash 

0.26 

0.67 

0.08 

0.20 

0.79 

0.71 


PPT: Precipitation. 

GOD: Growing degree days. 

PET: Potential evapotranspiration based on free water surface evaporative demand. 
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Temperature (°C) 

Fig. 6. Thermal gravimetric and derivative thermal gravimetric analysis of big 
bluestem. 

of big bluestem ranged from 17.28 to 19.05 MJ/kg. Prediction model 
1 (Dulong-Bertholot equation) overpredicts the HHV of samples 
and appears to have bigger bias error to experimental value. 
Nevertheless, the deviation of the prediction model 2 based on 
simple and rapid proximate analysis falls in a narrower range with 
the least RSS (residual sum of squares), indicating proximate 
analysis is good choice in terms of the accuracy of predicting HHV 
(Table 4). According to the correlations between HHV and proxi¬ 
mate analysis components (volatile matter, fixed carbon, and ash) 
in this study, the resulting Equation 5 has been derived from 
multiple linear regression analysis with factors (volatile matter and 
fixed carbon) using a least squares fitting program. This equation 
accounted for 85.5% of variation with a standard error of the esti¬ 
mate of 254.68: 

High heating value = -1109.78 + 207.46 volatile matter 

+ 210.82 fixed carbon (5) 

In general, higher fixed carbon and volatile matter contents 
resulted in a higher HHV. Big bluestem from Illinois had signifi¬ 
cantly higher HHV (average of 18.67 MJ/kg) than those from Kansas 
(average of 17.64 MJ/kg for Colby, KS; average of 17.73 MJ/kg for 
Hays, KS; and average of 18.07 MJ/kg for Manhattan, KS) (Table 4). 
In addition to location, the effect of ecotype on HHV showed the 
same trend as on specific heating of big bluestem. KAW had 
significantly higher HHV than other ecotypes, which might be 
partly attributed to the higher carbohydrate contents of native 
cultivar. In general, HHV of big bluestem was significantly affected 
by location (p < 0 . 01 ) and ecotype (p < 0 . 01 ), with the former being 
more influential (as shown by a large F-value). The interaction ef¬ 
fect of ecotype and location on HHV was statistically insignificant 
(p > 0.05). 

Table 4 shows the results of the proximate analysis calculated on 
a dry basis. The volatile matter content of big bluestem ranged from 
68.75% to 74.65%. The volatile matter was exclusive of moisture 
vapor and consists of permanent gases such as CH 4 , CO 2 , and CO 
forming the bio-oil after condensation [35]. Fixed carbon, produc¬ 
ing a char and burning as a solid material in the combustion system, 
ranged from 16.75% to 19.24%. Ash, the inorganic residues, ranged 
from 4.35% to 11.46%. Variations in the fixed carbon, volatile, and 
ash contents among the 16 samples were analyzed by two-way 
ANOVA to examine the genetic and environmental effects on 
proximate analysis of the big bluestem. With the exception of the 
interaction effect on the fixed carbon, significant effects of planting 
location, ecotype, and an interaction between location and ecotype 


were observed on the volatile matter, fixed carbon, and ash con¬ 
tents. Location had larger F values (50.16-1426.49) than ecotype 
(5.51-278.24) and interactions (1.53-85.55), showing that location 
effects were always highly significant, with much larger F values on 
proximate analysis of big bluestem, an order of magnitude larger 
than F values of ecotype and their interaction for volatile matter 
and fixed carbon contents, and approaching more than three orders 
of magnitude larger for ash content (Table 5). Table 5 shows envi¬ 
ronmental factors, including mean annual precipitation since 1961, 
potential evapotranspiration, and aridity index variation, which 
had a significant effect on the heat value with large coefficients of 
determination (F 2 >0.90). The mean annual precipitation since 
1961 in Illinois is twice as high as that in Colby, which provides a 
better environment for big bluestem with high heating value. For 
proximate analysis, variables were also associated with environ¬ 
mental predictors; the potential evapotranspiration explained 64%, 
86 %, and 79% of the variation in volatile matter content, fixed car¬ 
bon content, and ash content, respectively. It is suggested that the 
potential evapotranspiration played a greater role in divergence of 
the high heating value and proximate analysis variables of the big 
bluestem. 

4. Conclusions 

Specific heat of big bluestem was significantly affected by 
planting location, ecotype, and interaction between location and 
ecotype, but location was the dominant factor. A positive correla¬ 
tion model was developed between specific heat and temperature. 
There was no significant effect on thermal conductivity of big 
bluestem. A positive relationship between thermal conductivity 
and density and a negative relationship between the particle size 
and thermal conductivity were observed. With the exception of 
weight loss, planting location alone had a significant effect on 
thermogravimetric parameters of big bluestem. As a general 
conclusion, both planting location and ecotype significantly 
affected high heating value, but the former was more influential. 
Potential evapotranspiration was the most significant environ¬ 
mental factor affecting all thermal properties of big bluestem 
selected for this study. 
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List of nomenclature 


C p : specific heat, kJ/kg/K 
Cps : specific heat of sapphire, kJ/kg/K 
E: specific heat calibration constant 
E re f voltage of the reference resistor, V 
H: measured heating value, mW 
H r : heating rate, K/min 
k: thermal conductivity, W/m/K 
M: sapphire mass, mg 

q: heater power dissipated per unit length, W/m 

Rm: resistance of probe heater, Q 

R re f .; resistance of the reference resistor, Q 

T: final temperature, K 

t: final time, s 

To: initial temperature, K 

to: initial time, s 

T en( i: final temperature of decomposition peak in DTG curve, °C 
Tmax'- peak temperature of decomposition peak in DTG curve, °C 
Tone- initial temperature of decomposition peak in DTG curve, °C 



